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! Abstract 

o . 

, The independent solutions of the one-dimensional Schrodinger equation are approximated 

CsJ ■ by means of the explicit summation of the leading constituent WKB series. The continuous 

. matching of the particular solutions gives the uniformly valid analytical approximation to the 

wave functions. A detailed numerical verification of the proposed approximation is performed 
, for some exactly solvable problems arising from different kinds of potentials. 
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1 Introduction 

Perturbation theory, the variational method and the WKB approximation are very extensively 
used in quantum mechanics. If we deal with perturbation theory or with the variational method 
then similar questions arise. How to find the unperturbed Hamiltonian or how to find the trial 
function for an arbitrarily given potential? Universal answers are absent. In this sense both men- 
tioned methods are incomplete. In contrast, the WKB approximation is directly determined by 
Qh a given potential. However the conventional WKB approximation has unphysical singularities. 

An old problem in semiclassical analysis is the development of global uniform approximations 
to the wave functions. In previous works [1, 2], an essential improvement of the WKB approach 



was introduced for the logarithmic derivatives of the wave functions. In the present paper, we 
construct the second-order continuous approximation to the wave functions. The quality of the 
approximate wave functions is verified by means of a comparison with the exact solutions for 
different kinds of potentials. 
| We consider the linear one-dimensional Schrodinger equation 

^ = fW>, a) 

where Q(q) = 2m(V(q) — E) for an arbitrary potential V(q). The logarithmic derivative 

dln^(q,h) 
Y ^ %)= dq 

of a wave function ^(q, Ti) satisfies the nonlinear Riccati equation 

dY(q,h) 2 Q(q) 

—^ + {Y{qM =^-. (2) 

The WKB approach deals just with functions Y(q, K). In this approach, two independent solu- 
tions Y^(q,H) of the Riccati equation are represented by their asymptotic expansions 

Y±(q,h) = h- 1 (±Q 1 / 2 + f] h n Y±(q)) (3) 
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in powers of Plank's constant fi. The usual WKB approximation contains a finite number of 
leading terms Y^(q) from the complete expansions Y^(q,h). This approximation is not valid 
at turning points where Q(q) = 0. 

As it is well known, the WKB series is divergent. Numerous references regarding asymptotic 
expansions may be found in [3]. The direct summation of a divergent series does not exist. 
By summing one means finding a function to which this series is the asymptotic expansion [4] . 
In recent years many studies have been devoted to extracting some useful information about 
the exact eigenfunctions from the divergent WKB series (see, e.g., [5] and references therein). 
There are several investigations on the properties of the WKB terms [6, 7]. Unlike entirely 
exact but very complicated methods for some classes of potentials (see, e.g., [8]) our new way 
of using the WKB series gives an approximate but very simple and universal method of solving 
the Schrodinger equation. 



2 Explicit summation of the constituent WKB series 

Since [2] is likely to be inaccessible for the large majority of readers we reproduce previous 
results. First of all, the analysis of the well-known recursion relations [4, 6] 

Y± +1 = -(2y ± )~ 1 (^2 Yt Y n + i-> + dY±/dq^J , Y± = ±Q^ 2 (4) 
shows that the WKB terms are of the form 

n 

Y^( q ) = Q^J2 A UQ'(^Q , '(^---Q {j) ('i))Q i ~ 1 ^ (5) 

3=1 

where Q'(q) = dQ{q)/dq, Q"{q) = d 2 Q{q)/dq 2 and Q^\q) = d^Q{q)/dq^. Second, the substitu- 
tion of (5) into (3) allows us to reconstruct the asymptotic WKB series as an infinite sum 

oo 

Y±(q,h) = ±h- 1 Q^+Y / Z± J (q,h) 

3=1 

of new constituent (partial) asymptotic series 



3n+l 



Zts,M,K) = {h^y- 2 Y. {-§^)~ 2 A± J (Q>(q),Q"(q),...Q^(q)) 



n=j 



m powers of the ratio Q/h 2/3 . With the help of the recursion relations (4) we derive simple 
expressions 

<i = (QTB^, < 2 = Q"(QT~ 2 B± 2 

for two leading sequences of coefficients A^j . Here the numbers are determined by the 
following recursion relations: 

and B^ 2 is connected with B^ tl as follows: 

B% 2 = -\nBi^ n>2. (7) 
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The complete series Y^(q, h) are approximated by a finite number of leading constituent se- 
ries Z^ s j (q, h) in contrast to the use of a finite number of leading terms (q) in the conventional 



WKB approach. If we can find functions Zf(q,Ti) which are represented by asymptotic expan- 



j 



sions Z^ s j(q,h), then we obtain new approximations to the solutions of the Riccati equation. 
The number of used constituent series corresponds to the order of a proposed approximation. For 
instance the expressions ±fi~ l Q 1 / 2 + Zf{q, ft) are interpreted as the first-order approximations. 
In this paper, we consider only the second-order approximations ±h~ 1 Q 1 / 2 + Z^ z (q, h) + Z^(q, fi). 
Introducing the dimensionless variable 

< %\ 1 g(g) (*\ 
we are able to rewrite the leading constituent expansions in the form 

ifr-V 72 + Zt tl (q,h) + Z^ 2 (q,h) = ^^-^(a) + |^ 2 (a), (9) 
where we separate the asymptotic series in a 

CO 

^>) = ±« 1/2 + E^V~^ (io) 

71=1 



71=2 



The leading terms 



±a i/2 _ L-\ (12) 

± i»" 3/2+ i»" 3 < 13 > 

of these series may be deduced by using equations (6) and (7). 

Our aim is to sum constituent series (10) - (11). In other words, we must find functions 
yf{a) which are represented by these expansions. In order to perform the identification we 
substitute the approximate function 

*i(**) = ^w^ yt{a) + % vt{a) (14) 

into the Riccati equation (2). As a result we get the following equations: 

d 4 + («t? = «. (15) 

d yf , o ± ± 1 ( n d yt ±\ / 1r n 

-fo + 2yi 2/2 = 3 r^r - yi j (i6) 

for the functions yf(a). Direct verification shows that the asymptotic expansions (10) and (11) 
satisfy these equations. 

Equation (15) is the Riccati equation for the logarithmic derivatives of linear combinations 
of the well-studied Airy functions Ai(a) and Bi(a) [9]. We select particular solutions by means 
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of the known asymptotics (12). In the classically allowed region where Q(q) < (a < 0) we 
derive the explicit expressions 



I/f(a) = — ln(Bi(a)=FiAi(a)) (17) 
and in the classically forbidden region where Q(q) > (a > 0) we get the other solutions 

yr(a) = ^ In Ai(a), y+{a) = ± lnBi(a). (18) 



Finally, we can obtain the solutions of the linear equation (16) with asymptotics (13) in the 
closed form 

..±f„\ _ 1 f o„2,..±r„\\2 a„„.±(„\ , o„3 



J/2 (a) = ^ [-8a 2 (yf (a)) 2 - Aayf (a) + 8a J - 3 
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(19) 
(20) 



Although the functions (17) - (20) have the asymptotic expansions (10) - (11) if \a\ is large 
it should be stressed that the obtained functions possess different expansions if \a\ is small. 
Replacing yf by yf in expression (14) we get the second pair 



Q' 



/j2/3 | Q, 1 2/3 



n" 



(21) 



of approximate solutions. 

It is not surprising that the asymptotics of our approximation coincide with the WKB 
asymptotics far away from the turning points. At the same time our approximation reproduces 
the known [4] satisfactory approximation near the turning points. Naturally, our approximation 
gives the exact result for the linear potential V(q) = kq of a uniform field. Note that this 
potential represents an example of the explicit summation of the WKB series for the logarithmic 
derivative of a wave function. 



3 Approximate wave functions for the two-turning-point prob- 
lem 

With the aid the uniformly valid approximation to solutions of the Riccati equation derived in 
preceding section, we can now construct approximate wave functions. We consider the problem 
with two real turning points g_ and q + (q + > q~). The potential has its minimum at point q m . 
The first and second derivatives of the smooth potential are continuous at point q m . 

Two pairs of independent solutions of the Schrodinger equation are approximated by func- 
tions 

*±(g)=exp (^J 9 Y±(q')dq'y 
=exp (J q Y±(q')dq'y 

In accordance with the requirements of quantum mechanics, we must retain only the decreasing 
solutions *a p (g) in the classically forbidden regions (q < g_ and q > q + ). In the classically 
allowed region (<?_ < q < q + ) we retain a linear combination of two oscillatory solutions *i,(g) 
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and ^ ap (q). By matching particular solutions at the turning points g_ and q + , we obtain the 
continuous approximate wave function which is represented by following formulas: 

= Ccos | exp (- j*~ Y~^)d/j (22) 

if q < qr_, 

Ml) = Cexp ^ dq j cos — dq - - j (23) 

if q- < q < q + , and 

*s(?) = C(-l)« cos I exp y-(^ exp (jT ^l^dg') (24) 
if g > q+. 

Here we have the new quantization condition 

f q+ Q' Y+(q,E)-Y-(q,E) j 2 

X_ 2i dq = ,(n+-), n = 0,1,2,... (25) 

which determines the spectral value E sp (n) of energy implicitly. We denote the wave functions 
with E = E sp (n) as ^> ap (q,n). Then we may choose the value of an arbitrary constant C in 
order to ensure the usual normalization ap (n)\^ ap (n)) = 1 where |\P ap (n)) is the vector in 
Hilbert space which corresponds to the function ^f ap (q,n). The proposed approximation is an 
alternative to the well-known [4, 10] Langer approximation [11] which employs a /i-expansion 
different from the WKB series. 

Thus the approximate eigenfunctions are determined completely. However a question arises 
regarding the optimal approximate eigenvalues, because the value E sp (n) is not a unique choice. 

Since explicit expressions for wave functions have already been obtained, we are able to 
calculate the expectation values 

E(n) = {* ap (n)\H\V ap (n)) (26) 

of the Hamiltonian 

In accordance with the eigenvalue problem 

H\V) - E\V) = 

we construct the discrepancy vector 

\D(e,n)) = H\^ ap {n)) - e, |* op (n)), 

where e is an arbitrary parameter while H and \^f ap (n)) are given. It is natural to require 
that the discrepancy vector should not contain a component proportional to the approximate 
eigenvector. In other words, we consider the orthogonality condition 

(y ap (n)\D(e,n))=E(n)-e = 

as a criterion for the selection of the optimal approximate eigenvalue. As a result we just get 
E(n) while E sp (n) does not fulfil the above requirement. It should also be noted that the scalar 
product (D(e, n)\D(e, n)) is minimized at e = E(n). 
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4 Verification of the proposed approximation 

Now we must verify our approximation numerically for exactly solvable problems. We compare 
the normalized approximate wave functions ^> ap (q,n) with the normalized exact wave functions 

In order to estimate the closeness of two functions fi(q) and /2(g), we consider two corre- 
sponding vectors and I/2) in Hilbert space. Then we construct a deviation vector |A/) = 
— I/2} and a scalar product 

(A/|A/) = {h\h) + </ 2 |/ 2 ) - {h\f 2 ) - {f 2 \h). 

Now we can define the relative deviation 

„_ 1 (/ll/2) + </ 2 |/l) 
/_ (/l|/l) + (/2|/ 2 ) (27) 

as a numerical estimate of the closeness of two functions. Note that 5f = if fi(q) = / 2 (<z)- 
Thus, we get the following estimate 

6®( n ) = 1 - <*ez(n)|* ap (n)) (28) 

in the case of the normalized real functions ^f ex (q,n) and ^ ap (q,n). The same numerical com- 
parison may be performed 

2(tf' (n)ltf' (n)) 

S^i'(n) = 1 p (29) 

{ ' (*tx(n)l* / exW> + <*^(n)l*^(n)> 1 ' 

for the first derivatives = d^>(q)/dq . Naturally, we can define analogous estimates for 

higher derivatives. 

Moreover, we estimate the closeness of functions H^f ap (q,n) and H^ ex (q,n) with the help 
of the following quantity 

sA y, n) = ^a P {n)\H 2 \^a P {n)) + E ex {n?{l - 2(* ex (n)\* ap (n))) 

^ ap (n)\H^ ap (n))+E ex (nf 

where E ex (n) is the exact energy value. 

In addition, we compare two functions H^ ap (q,n) and E(n)^f ap (q,n). This comparison 
may be performed when we do not know the exact solutions. As a result we get the relative 
discrepancy 

d(n) = ^a P (n)\H 2 \^ ap (n))-E(nf 
{m ap {n)\H^ ap (n)) + E{nf 

which is directly connected with the Schrodinger equation under consideration. 
Finally, we characterize our approximation by the usual relative energy error 

The verification is performed for three potentials with different asymptotics. They are the 
harmonic oscillator potential 

V(q) = kq 2 , 
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the Morse potential 

^ («-*'' -*-) 

and the modified Poschl- Teller potential 

= A(A-l)fiV 
2m cosh 2 (ag) 

The exact solutions of the Schrodinger equation for these potentials may be found in [12]. 

Table 4.1 Numerical verification of the proposed approximation. 
n 5^(n) W(n) 5HV(n) d{n) AE(n) 



The harmonic oscillator potential. 






2.16 


io~ 4 


4.60 • 10~ 3 5.79 • 10~ 2 


5.32 


10" 2 


4.52 


10" 


-3 


1 


1.59 


10~ 5 


7.12 • 10~ 5 1.33 • 10" 4 


6.09 


IO" 5 


5.64 


10" 


-5 


2 


3.86 


10~ 6 


1.05 • 10~ 5 1.48 • 10~ 5 


3.78 


10~ 6 


7.15 


10" 


-6 


3 


1.50 


IO" 6 


3.10 • 10~ 6 4.12 • IO" 6 


7.31 


io- 7 


1.89 


10" 


-6 








The Morse potential 


( 7 = 4.5). 











2.76 


10" 4 


5.98 • 10~ 3 4.59 • 10" 3 


5.77 


10~ 3 


-1.45 


io- 


-s 


1 


4.95 


10~ 5 


2.34 • 10" 4 2.45 • 10" 4 


3.56 


10~ 4 


-1.60 


io- 


-4 


2 


2.26 


10~ 5 


8.22 • 10~ 5 3.72 • 10~ 4 


4.51 


10~ 4 


-1.01 


io- 


-4 


3 


1.77 


10~ 5 


6.43 • 10~ 5 2.35 • 10~ 3 


2.49 


IO" 3 


-1.59 


10" 


-4 






The modified Poschl- Teller potential 


(A = 


5). 









2.31 


10" 4 


6.60 • 10~ 3 4.49 • 10~ 3 


5.70 


10" 3 


-1.44 


10" 


-3 


1 


2.91 


10~ 5 


2.10 • 10" 4 1.65 • 10" 4 


2.97 


io~ 4 


-1.62 


io- 


-4 


2 


1.98 


10~ 5 


6.17 -10~ 5 2.36 -10" 4 


3.18 


10" 4 


-1.01 


io- 


-4 


3 


5.22 


lO" 5 


4.90 • IO" 5 1.95 • IO" 3 


2.05 


IO -3 


-1.48 


10" 


-4 



Table 4.1 shows that our approximation gives fairly accurate results for all considered po- 
tentials and for all considered quantities. Hence, the performed reconstruction of the WKB 
series and subsequent explicit summation of the leading constituent (partial) series yield the 
satisfactory (qualitative and quantitative) description of wave functions. 
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